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We give two combinatorial interpretations of the Matrix Ansatz of the PASEP in terms of lattice paths and rook 
placements. This gives two (mostly) combinatorial proofs of a new enumeration formula for the partition function of 
the PASEP. Besides other interpretations, this formula gives the generating function for permutations of a given size 
with respect to the number of ascents and occurrences of the pattern 13-2, the generating function according to weak 
exceedances and crossings, and the n* moment of certain g-Laguerre polynomials. 

Resume. 

Nous donnons deux interpretations combinatoires du Matrix Ansatz du PASEP en termes de chemins et de placements 
de tours. Cela donne deux preuves (presque) combinatoires d'une nouvelle formule pour la fonction de partition du 
PASEP. Cette formule donne aussi par exemple la fonction generatrice des permutations de taille donnee par rapport 
au nombre de montees et d'occurrences du motif 13-2, la fonction generatrice par rapport au nombre d'excedences 
faibles et de croisements, et le n'™'^ moment de certains polynomes de g-Laguerre. 
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1 Introduction 

In recent work of Postnikov [p^, permutations were given a new description as pattern-avoiding fill- 
ings of Young diagrams. More precisely, Postnikov made a correspondence between positive Grassmann 
cells, pattern-avoiding fillings called J-diagrams, and decorated permutations (which are permutations 
where the fixed points are bi-coloured). In particular, the usual permutations are in bijection with per- 
mutation tableaux, a subclass of J-diagrams. Permutation tableaux have subsequently been studied by 
Steingrimsson, Williams, Burstein, Corteel, Nadeau [Q, 0, ||, |8]|, and proved to be very useful for work- 
ing on permutations. 
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Rather surprisingly, Corteel and Williams established a link between permutation tableaux and the sta- 
tionary distribution of a classical process studied in statistical physics, the Partially Asymmetric Exclusion 
Process (PASEP). This process is described in ||^, Briefly, the stationary probability of a given state 
in the process is proportional to the sum of weights of permutation tableaux of a given shape. The factor 
behind this proportionality is the partition function, which is the sum of weights of permutation tableaux 
of a given half-perimeter 

An alternative way of finding the stationary distribution of the PASEP is given by the Matrix Ansatz 
Suppose that we have operators D and E, a row vector {W\ and a column vector \ V) such that: 

DE-qED = D + E, {W\E = {W\, D\V)^\V), and {W\\V) ^ I. (1) 

Then, coding any state of the process by a word w of length nin D and E, the probability of the state w 
is given by (W^juilF) normalised by the partition function + V^). 

We briefly describe how the Matrix Ansatz is related to permutation tableaux [^. First, notice that there 
are unique polynomials Uij e Z[g] such that 

(D + E)" = J2 n^,3E'D^ 

i,j>0 

This sum is called the normal form of (_D+i?)". It is useful since, for example, the sum of coefficients Uij 
gives an evaluation of (T/F|(-D + E)'"-\V). Each coefficient Uij is a generating function for permutation 
tableaux satisfying certain conditions, or equivalently, alternative tableaux as defined by Viennot [^]. 

We give here two combinatorial interpretations of the Matrix Ansatz in terms in lattice paths and rook 
placements, and get two semi-combinatorial proofs of the following theorem: 

Theorem 1 For any n > 0, we have: 



{W\{yD+Ey^-^\V) = ^^(-1)'= - E 



fc=0 \ j=0 




The combinatorial interpretation of this polynomial, in terms of permutations, is given in Proposition 
For y — I this specialises to: 

Corollary 1 For any n > 0, we have: 

{W\{D + Er-^\V) = -—-Y,[-lf[{^ 



2n \ 
i-k) 




Besides the references mentioned earlier, we have to point out an article of Williams |]27|], where we 
find the following formula for the coefficient of y'"^^ in {W\{yD + E)'^\V): 

EmAl) = E - *]?'?™""' (C)?""^ + (.-l)) ■ (2) 



It was obtained by enumerating J-diagrams of a given shape and then computing the sum of all possible 
shapes. Until now it was the only known polynomial formula for the distribution of a permutation pattern 
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of length greater than two (See Proposition p]). A lthough the article focuses on J-diagrams, Williams 
and her coauthors sketched in Section 4 of [|14{] how this could have been done directly on permutation 
tableaux. Recently, Williams's formula has been obtained also by Kasraoui, Stanton and Zeng in their 
work on orthogonal polynomials [|l2|. We will show in the last section how our formula can be applied to 
prove and extend a conjecture presented in [p7|]. 

The polynomial y{W\{yD + E)'"'^^\V) was already heavily studied. 

Proposition 1 For any n> \ the following polynomials are equal: 

• y{W\{vD + EY-^\V), 

• the generating function for permutation tableaux of size n, the number of lines counted by y and 
the number of superfluous 1 's counted by q |26|/, 



• the generating function for permutations of size n, the number of ascents counted by y and the 
number of 13-2 patterns counted by q ^ [7i| /, 

• the generating function for permutations of size n, the number of weak exceedances counted by y 
and the number of crossings counted by q ^ 



• the generating function ofPDSAWs ( partially directed self-avoiding walks ) in the asymmetric wedge 
of length n where the number of descents is counted by y and the number of north steps is counted 
byqfy. 



• the n moment of the q-Laguerre polynomials / |72[ |27| /. 

Remark. We can view the formula in Corollary |l] as an analog of the Touchard-Riordan formula for 
the number of matchings of 2n according to the number of crossings: 



M matching of 2n \ i / j,^q \ \ / \ 



q 2 . 



We remark that this formula also gives the 2n* moment of the q-Hermite polynomials. 

In pcj], Penaud gave a combinatorial proof of this formula. By generalising Penaud's method we 
conjectured Theorem ^ and were hoping for a completely combinatorial proof thereof. However, at the 
time of writing the last step of this combinatorial proof is still missing. 

This article is organised as follows: we first show how the Matrix Ansatz is naturally related to lattice 
paths. Then we give two proofs of our main Theorem, one based on lattice paths and the other one based 
on rook placements. We end with a discussion and some applications. 

2 A first proof using lattice paths and functional equations 

2. 1 The Matrix Ansatz and lattice paths 

We follow the ideas developed in ||]. Looking for a solution of the system defined in Equation (|l]) we 
find: 
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and E = 








.. + q' 


^hi = 


[I- 


.. + q' 


{W\ = 


(1,0,0, 


. . .), and 


1^) = 


(1,0,0, 


...f . 



>o such that 

ifi equals j — 1 or j, 
otherwise, 

ifi equals j or j + 1, 
otherwise, 



Then these matrices and vectors satisfy the Ansatz of Equation 

We can interpret y{W\{yD + E)"^^\V) as the generating polynomial of paths of length n — 1. The 
weight of a path is the product of the weight of its steps and the weight of the starting and ending points. 
If a path starts (resp. ends) at (0, i) (resp. {n — 1, i)) the weight of the starting (resp. ending) point is Wi 
(resp. Vi). The weight of a step going from {x, i) to {x + l,j) is Dij + Eij. We call i the starting height 
of the step. See ^ for details. 

Proposition^ implies that the paths we are dealing with here are bi-coloured Motzkin paths, i.e., paths 
that start and end at height zero and consist of north-east, south-east and two types of east steps. Using a 
classical bijection we can transform these paths of length n — 1 into Motzkin paths of length n where east 
steps of type 2 can not appear at height zero. 

Propositions y{W\{yD + E)"'^^\V) is the generating polynomial of weighted bi-coloured Motzkin 
paths of length n such that the weight of steps starting at height i is 

• y + yZ + • • ■ + VQ^ = y ^i-q fa'' north-east steps and east steps of type 1, and 

• 1 + g + . . . + q^^^ = for south-east steps and east steps of type 2. 
This can also be done combining results in |[ p^ . 

2.2 The proof 

The method used in this subsection is inspired by an article of Penaud [^o|]. We extract a factor of (1 — g)" 
from the generating polynomial of the weighted bi-coloured Motzkin paths from Proposition |] and obtain 
that 

y{W\{yD + Er-^\V) ^ J2 Hp), 

where P{n) is the set of labelled bi-coloured Motzkin paths of length n such that the weight of steps 
starting at height i is either 

• y or —yq^^^ for north-east steps or east steps of type 1. 

• 1 or — for south-east steps or east steps of type 2. 

Let M{n) be the subset of the paths in P{n) such that the weight of any east step and the weight of any 
peak (a north-east step followed by a south-east step) is neither 1 nor y. Let Ain,k,j be the number of left 
factors of bi-coloured Motzkin paths of length n, final height k, and with j south-east and east steps of 
type 1. 
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Lemma 1 There is a bijection between paths in P{n) and pairs of paths such that 

• the first path is a left factor of a bi-coloured Motzkin path of length n and final height k 

• the second path is in M{k), for some k between and n. 
In particular, we have 

n n — k 

peP{n) fe=Oj=0 p£M{k) 

Proof: Let p be a path in P{n). We decompose p into a sequence mi gi 771252 • • • mkqkfnk+i such that 

• the nii are maximal (but possibly empty) sub-paths of p with all steps having weight 1 or y, and 
returning to their starting height, 

• the qi are single steps. 

It follows that q\q2 . . . gfe is a path in M{k). Replacing in the sequence migi 772292 • • • nikqkmk+i each 
step qi by a north-east step, and taking into account the number of south-east steps and east steps of type 1 , 
we obtain a path in A4n,k,j of weight . □ 

It remains to compute Mn,k,j and Mk = J2peM(k) ^(^)- 
Proposition 4 The number M.n,k,j of left factors of bi-coloured Motzkin paths of length n, final height 
k, and with j south-east steps and east steps of type 1, is (^) {^^^ — (j+fc+i)- 

Proof: We note that the formula can be seen as a 2 x 2 determinant. By the Lindstrom-Gessel-Viennot 
lemma, this equals the number of pairs of non-intersecting lattice paths taking north and east steps from 
(1, 0) to (ti — j, j) and (0, 1) to {n — j — k, j -\- k) respectively. 

We transform such a pair of paths step by step into a single Motzkin path according to the following 
translation table: 



step of 


lower path 


upper path 


Motzkin path 




north 


north 


east type 1 




east 


east 


east type 2 




north 


east 


north-east 




east 


north 


south-east. 



It is easy to see that the condition that the two lattice paths do not intersect corresponds to the condition 
that the Motzkin path does not run below the x-axis. Furthermore, we see that the number of east and 
south-east steps equals j, the number of north steps of the lower path. □ 

Proposition 5 The generating polynomial M). equals Yli=o y'q'*^*'"'"^"'^. 

Proof: We add an extra parameter on the paths in M{n), that marks the number of steps that have a 
weight different from 1 and y. More precisely, the weight of steps starting at height i is 

• y or —yzq^~^^ for north-east steps or east steps of type 1, and 
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• 1 or — zg* for south-east steps or east steps of type 2. 

Let M (z) = X)n>o ^" X)peA/(n) ^(p)- obtain a functional equation for M{z) by considering 

the following decomposition: A path is either (a) empty, (b) a north-east step of weight 1, followed by 
a path, followed by a south east step of weight —qyz, followed by another path, (c) a north-east step of 
weight —qz, followed by a path, followed by a south east step of weight —qyz, followed by another path, 
(d) a north-east step of weight —qz, followed by a path, followed by a south east step of weight y, followed 
by another path, (e) a north-east step of weight 1, followed by a non-empty path, followed by a south east 
step of weight y, followed by another path, (f) an east step of type 1 followed by another path, or (g) 
an east step of type 2 followed by a path. The corresponding weight is (a) 1, (b) —M{qz)qyzM{z)t^, 



(c) qzM{qz)qyzM{z)t^, (d) -qzM{qz)yM{z)t^, (e) {M{qz) 
—zM{z)t, respectively. Thus, we have: 

M{z) = 1- {qyzt + zt + yt^)M{z) + yt'^{l- 
PTOcseding similar to [p^, we use the linearising Ansatz 

1 H{qz) 



1) yM{z)t\ (f) -qyzM{z)t, or (g) 
qzfM{z)M{qz) . 



M{z) 



1 - z H{z) 



to obtain 



H{z) - (1 + yt^)H{qz) + yt^H{q^z) = z {H{z) + (1 + qy)tH{qz) + qyt^H{q^z)) . 

Solving recursively for the coefficients c„ of H{z) = J2'^=o ^nz", we obtain a solution in terms of a basic 
hypergeometric series. 



H{z) = 201 (-i, -tqy; t'^qy; q, z) 



1^— 



Note that we are dealing with a series of the type 20i (a, b; ab; q, z) where a 
to take the limit z ^ 1, we need to transform using Heine's transformation 



-t and h = —tqy. In order 



We find that 



and therefore 



2<^i(a, b, ab; g, z 
M{z) = 



(az, 6; q)r 



24>i{a,z;az;q,b) . 



{ab,z;q)c 

1 2(f>i{a,qz;aqz;q,b) 
1 ~ az 2</'i(a, az; g, &) 



M(l) 



1 



1 - a 

Changing back to a = —t and b = —tqy. 



>0i(a,g; aq;q,b) 



n=0 



aq" 



m-\-n—k i—{ 



i{k-i+l) 



Combining the previous results, we get a proof of Theorem [TJ 



□ 
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3 A second proof using the Matrix Ansatz and rook placements 



For further details about material in this section, see Ql One of the ideas at the origin of this proof is 
the following. From D and E of the Matrix Ansatz, we define new operators D and E as 

(7-1 1 - - 1 1 

D = - D+-, E = - E+-. 

q q q q 

An immediate consequence is that 

DE^qED = ^-^, (WlE^lWl, and D\V) = \V). (3) 

q 

This commutation relation is somewhat simpler than the one satisfied by D and E, as it has no terms 
linear in D or E. Moreover, we have q{yD + E) + {1 — q) {yD + E) = l + y, for any parameter y. Using 
this identity, we obtain the following inversion formulae between {yD + i?)" and {yD + E)": 

{^-qr{yD + Er = J2(''^i^ + yr''i-^)'q'iyD + E)^ and (4) 

k=0 ^ ^ 



q"{yD + ET = E \k) + - + Ef- (5) 

In particular, the first formula means that in order to compute the coefficients of the normal form of 
{yD + E')", it is sufficient to compute the ones of {yD + E)^ for all < fc < n (as taking the normal 
form is a linear operation). 

Except for a q-dependent factor, the operators D and E are also defined in [ ^3| ] and In the first ref- 
erence, Uchiyama, Sasamoto and Wadati used the commutation relation between D and E to find explicit 
matrices for these operators. They derive the eigenvalues and eigenvectors of D + E, and consequently the 
ones of D + E, in terms of orthogonal polynomials. In the second reference, Blythe, Evans, Colaiori and 
Essler also use these eigenvalues and obtain an integral form for (£> + E)"\V). They also provide an 
exact integral-free formula of this quantity, although quite complicated since it contains three summations 
and several g-binomial coefficients. 

In this article, instead of working on representations of D and E and their eigenvalues, we study the 
combinatorics of the rewriting in the normal form of {D + i?)", and more generally {yD + E)" for 
some parameter y. In the case of D and E, the objects that appear are the rook placements in Young 
diagrams, long-known by combinatorists since the results of Kaplansky, Riordan, Goldman, Foata and 
Schiitzenberger (see [ p^ and references therein). This method is described in [^, and is the same that 
the one leading to permutation tableaux or alternative tableaux in the case of D and E. 

Definition 1 Let X be a Young diagram. A rook placement of shape A is a partial filling of the cells of X 
with rooks ( denoted by a circle o ), such that there is at most one rook per row ( resp. per column ). 

For convenience, we distinguish with a cross ( x ) each cell of the Young diagram that is not below (in the 
same column) or to the left (in the same row) of a rook (we are using the French convention). The number 
of crosses is an important statistic on rook placements, which was introduced in [|o| as a generalisation of 
the inversion number for permutations. Indeed, if A is a square of side length n, a rook placement R with 
n rooks may be visualised as the graph of a permutation a E 6„, and in this interpretation the number of 
crosses in R is the inversion number of a. 
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Definition 2 The weight of a rook placement R with r rooks, s crosses and t columns is w{R) ~ p^q'^y*, 
where p = 

With the definition of rook placements and their weights we can give the combinatorial interpretation 

of {W\{yD + E)"\V). 

Proposition 6 For any n, {W\{yD + is equal to the sum of weights of all rook placements of 

half-perimeter n. 

The enumeration of rook placements leads to an evaluation of {W\{yD^E)'^~^\V) , hence of {W\{yD+ 
E)^~^\V) via the inversion formula (Q). 

3. 1 Rook placements and involutions 

Given a rook placement R of half-perimeter n, we define an involution a{R) by the following construc- 
tion: label the north-east boundary of R with integers from 1 to n. This implies that each column or row 
has a label between 1 and n. If a column, or row, is labelled by i and does not contain a rook, it is a fixed 
point of a{R). Also, if there is a rook at the intersection of column i and row j, then a{R) sends i to j 
(and j to i). 

Given a rook placement R of half -perimeter n, we also define a Young diagram f3{R) by the following 
construction: if we remove all rows and columns of R containing a rook, the remaining cells form a Young 
diagram, which we denote by We also define ^(i?) = (a(i?), f3{R)). See Figure |l] for an example. 



R 



X 


X 








o 


X 


X 


X 










o 


X 






o 




X 




X 






X 



</>(i?) 



Fig. 1: Example of a rook placement and its image by the map (f>. 



Proposition 7 The map cf) is a bijection between rooks placements in Young diagrams of half -perimeter n, 
and ordered pairs {I, A) where I is an involution on {1, ... ,n} and A a Young diagram of half-perimeter 
\Fix{I)\. If <j){R) = (/, A), the number of crosses in R is the sum of \\\ and some parameter ^{I). 

Proof: This kind of bijection rather classical, see for instance |Q, |I3|]. Note that the pairs (/, A) may 
be seen as involutions on {1, . . . , n} with a weight 2 on each fixed point. For the second part of the 
proposition, we just have to distinguish different kinds of crosses in the rook placement R. For example, 
the crosses with no rook in the same line and column are enumerated by | A|. □ 

Corollary 2 Let Tj,k,n be the sum of weights of rook placements of half perimeter n, with k lines and j 
lines without rooks. Then for any j,k,n, we have: 

y^To^k-3,n- (6) 

Proof: The previous proposition means that the number of crosses is an additive parameter with respect 
to the decomposition R (/, A). This naturally lead to a factorisation of the generating function. □ 



j,k,n 



n-2k + 2j 
j 
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3.2 The recurrence 

Proposition 8 We have the following recurrence relation: 

To,k,n = To,fc,n-l +py[n+ 1 - 2k]qTo^k-l,n-l- (7) 

Proof: We have the relation: 

TQ.k,n = To,fc,n-l + pTl^k,n-l- (8) 

Indeed, we can distinguish two cases, whether a rook placement enumerated by To,k.n has a rook in its 
first column or not. These two cases give respectively the two terms of the previous identity. To end the 
proof we can apply identity to the second term. □ 

The recurrence is solved by the following formula. 

Proposition 9 

k 

(9) 



To,k.n = q^''J2(-^y'i'^ 



1=0 



n ~ 2k + i 
i 



n \ I n 

k — il \k — i — 1 



From this proposition, identity (||), and a q-binomial identity, we derive a formula for Tj ^ ^. 
Proposition 10 

k k 

Summing this identity over k gives the following result. 
Proposition 11 

{W\{yD + E)\V) = {1 + y)G{n) - Gin + 1), (10) 
where G(n) = E ( ( ■) " ! i) ) E .'^-^.^^'--^-^^ • 

This formula is a linear combination of the polynomials Pk = X]^=o q^^^^^^'^\ the coefficients being 
polynomials in y, just as in Theorem |l|. With this result and the inversion formula (^, we can prove 
Theorem[l|: the last step is a simple binomial simplification. 



4 Applications 

Among all the objects of the list in Proposition |l], the most studied are probably permutations and the 
pattern 13-2, see for example [|[ ^ |l8[ |l^. In particular, in |||, |l9|] we can find methods for obtaining, 
as a function of n for a given fc, the number of permutations of size n with exactly k occurrences of 
pattern 13-2. By taking the Taylor series of (|l]), we obtain direct and quick proofs for these results. As an 
illustration we give the formulae for A; < 3 in the following proposition. 
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Proposition 12 The order 3 Taylor series of {W\{D + Ey^^^\V) is: 

where C„ is the nth Catalan number. 

More generally, a computer algebra system can provide higher order terms, for example it takes no 
more than a few seconds to obtain the following closed formula for [q^^] {W\{D + 

io!(«+i2')!(n-8)! + ''^ + 2093 + 32354 7?° + 228543 - 318990 

-17493961 n^- 104051458 - 6828164 + 2022876520 71"* 
+6310831968 rt' + 5832578304 + 14397419520 n + 574801920o) , 



which is quite an improvement compared to the methods of [|19p. In addition to exact formula, we can give 
asymptotic estimates, for example for the number of permutations with a given number of occurrences of 
pattern 13-2. 



Theorem 2 For any fixed m > 0, 



[g"](W^|(Z? + ^)"-i|y) 



as n ^ 00. 



Proof: When n ^ 00, the numbers - („_^^_2) dominated by the Catalan number ;;;^(^")- 

This implies that in (1 — q)'"-{W\{D + Ey"-^^\V), each higher order term grows at most as fast as the 
constant term C„. On the other side, the coefficient of in (1 — g)^" is asymptotically n™ /m\. □ 

Since any occurrence of the pattern 13-2 in a permutation is also an occurrence of the pattern 1-3-2, a 
permutation with k occurrences of the pattern 1-3-2 has at most k occurrences of the pattern 13-2. So we 
get the following corollary. 

Corollary 3 Let ipkin) be the number of permutations in 6„ with at most k occurrences of the pattern 
1-3-2. For any constant C > 1 and fc > 0, we have 

*^<»> ^ W 

when n is sufficiently large. 

So far we have only used Corollary |l} Now we illustrate what can be done with the refined formula 
given in Theorem n] For example, when q = then the coefficient of ?/™ is given by the expres- 

sion ELo(-l)' C^DLlk) - This is equal to the Narayana number N{n,fn) = 

n (m) (m-i) ^^^^ ^ Combinatorial proof). 

We can also get the coefficients for higher powers of q. For example it is conjectured in p7[ ] that the 
coefficient of qy™ in {W\y{yD + i?)"^* |y) is equal to {^^i) („j'l2) ■ ^PPlyii^S '^^^ results we can prove: 
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Proposition 13 The coefficients ofqy"^ and q^y™" in {W\y{yD + EY~^\V) are respectively: 

/ ''I \( ''^ \ . + + nm + m — w? — 4 

\m+l)\rn-2) ^" \m-2)\m + 2) 2(n + l) ' 

Proof: A naive expansion of the Taylor series in q gives a lengthy formula, which is simplified easily after 
noticing that it is the product of (^) and a rational fraction of n and m. □ 
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